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Abstract 

The spectrum of primordial perturbations obtained by calculating 
the quantum gravitational corrections to the dynamics of scalar per¬ 
turbations is compared with Planck 2013 and BICEP2/iCec/: Array 
public data. The quantum gravitational effects are calculated in the 
context of a Wheeler-De Witt approach and have quite distinctive fea¬ 
tures. We constrain the free parameters of the theory by comparison 
with observations. 


1 Introduction 

The spectra of the inhomogeneities in the cosmic microwave background 
(CMB) contain much information on the evolution of the Universe from its 
early stages until now. The latest Planck mission [T] measured such inhomo¬ 
geneities and their statistical properties with great accuracy thus improving 
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the present constraints on cosmological evolution. Planck results indicate 
that the so called Standard Cosmological model can £t CMB data well on 
assuming that the spectra of the primordial cosmological perturbations at 
the onset of the hot big bang are given by "P ~ A{k/k^)'^ where is a 
suitable pivot scale. Such a power law dependence and a red tilted spectrum 
in the scalar sector may be produced by an inflationary era [2] driven by a 
single, minimally coupled scalar held with a hat potential |3j. 

Planck data conhrm the mild evidence, already hinted at by WMAP j3] 
and other independent datasets, of a power suppression of the temperature 
anisotropies spectrum on large scales. Despite that this evidence may well 
be explained as an ehect of cosmic variance, it is still worth investigating the 
intriguing possibility that it is associated with some physical ehect. 

Various attempts have been made in order to physically justify the power 
suppression on large scales [5]. In particular, it was observed that quantum 
gravitational ehects may ahect the large scales of the primordial spectra of 
perturbations [U]. Such an ehect is naively explained by the fact that the 
largest scales we observe in the CMB are those which hrst exit the horizon 
during inhation (and re-enter at later times) and are thus more ahected by 
large energy/curvature ehects than the shortest scales. 

Similar models for the evolution of primordial cosmological perturbations in 
the framework of quantum cosmology have been investigated recently jS]. In 
particular in the quantum gravitational corrections to the single, mini¬ 
mally coupled scalar held inhation are obtained through a Born-Oppenheimer 
decomposition [8j and without specifying any further assumption on the inha¬ 
tionary potential except for the requirement that it determines the slow-roll 
(SR) evolution. These corrections ahect the infrared part of the spectra and 
lead to an amplihcation or a suppression depending on the vacuum choice, 
the value of the SR parameters and the number of e-folds of inhation. More 
interestingly they depend on the wavenumber k and scale as k~^ in both the 
scalar and the tensor sectors. 

In this paper the predictions of this model are compared with observations 
and the free parameters of model are constrained. The paper is organized 
as follows: in Section 2 we review the basic formalism introduced in m ; in 
Section 3 we present the modihed primordial spectra obtained on calculating 
the quantum gravitational corrections; in Section 4 we discuss our analysis, 
we specify the dataset compared with the theoretical predictions and the free 
parameters involved in data fitting, we report and comment the Monte Carlo 
results; hnally the conclusions are drawn in Section 5. 
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2 Formalism 


The quantum gravitational effects on the spectra of the primordial perturba¬ 
tions Vk [9] have been obtained in a Wheeler-De Witt approach to quantum 
gravity ra. These effects can be evaluated to the first order in Mp 
solving the following differential equation 


-2 


on 


+ 2—p + = 0 

dr]'* dp dp 

for the two-point function 

Piv) = (0|f)fc|0) 

where the quantum gravitational corrections are given by 
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p is the conformal time, a is the scale factor and a prime means a derivative 
w.r.t. p. In the scalar sector Vk is the canonically quantized Mukhanov vari¬ 
able, = /c^ —^ with 0o is the homogeneous inflaton field and H 

is the Hubble parameter. In the tensor sector Vk is the canonically quantized 
rescaled tensor perturbation and uP = Let us note that the above 

master equation Q is only valid within a perturbative approach and Ap is 
the expression for the quantum gravitational corrections up to order Mp“^ 
to the evolution of the Bunch-Davies (BD) vacuum |11| . The quantum grav¬ 
itational correction term Ap is associated with “non-adiabatic” transitions 
between eigenstates of the invariant vacuum (BD) and higher eigenstates of 
the invariants m 

For realistic inflationary models, with a slow rolling inflaton field, the equa¬ 
tion ([^ can be further simplified to the first order in the SR approximation 
and one has 

2 (l -I- “iesR — Ipsr) 


u^ = k^- 

for the scalar perturbations and 
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for the tensor perturbations with 


esR = -^ and tjss 



( 6 ) 


where a dot means a derivative w.r.t. the proper time dt = adr]. Given the 
expressions 0.0 we observe that the results for the scalar sector can be 
transformed into the results for the tensor sector by simply replacing rjsR 
with esR. 


3 Quantum gravitational corrections 


The effect of Ap on the evolution of the two-point function p is that of adding 
to the standard BD solution prd (leading to the standard power-law expres¬ 
sion for the primordial spectra) a contribution of order Mp“^. Furthermore, 
to first order in the quantum gravitational corrections, the perturbative ap¬ 
proach adopted still leaves the freedom of considering a slightly modified BD 
prescription providing the standard vacuum is recovered when Mp —)• oo. In 
the scalar sector the final expression for p has the following form in the long 
wavelength limit: 


(L) (L) 

p'' > ~ ’ 
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where k is an unspecified reference wave number. The appearance of k in 
the quantum corrections can be traced back to the three volume integral 
in the original action for the homogeneous inflaton-gravity system plus per¬ 
turbations [7j. Such a volume, on a spatially flat homogeneous space-time, 
is formally infinite and consequently the value of k remains undetermined. 
Naively one may argue that k is related to an infrared problem (divergence) 
and indeed, in the literature, its value is taken to be the infrared cut-off 
for the perturbations, namely the largest observable scale in the CMB. The 
parameter 5k in ([^ is related to the initial conditions for the quantum state 
describing the evolution of the perturbations. Such initial conditions are 
fixed by the BD prescription at the unperturbed level but in principle a de¬ 
parture from pure BD of the order Mp“^ is consistent with the approach. In 
general 5k may depend on k and, in particular, on setting the SR parameters 
to zero one recovers the exact BD solution in de Sitter for 5k = —k/2. 
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The time dependent contribution —kt] in ([^ is evaluated at the end of in¬ 
flation. The wavenumber k necessarily refers to the scales, around the pivot 
scale /c*, which are probed by CMB and exited from the horizon iV* ~ 60 
e-folds before inflation ends. Its contribution to Q is 








OiendH j 


-“^{^SR-IISr) 


g2V*(esjj-?75fl) 
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and may well lead to a contribution of O {1) for reasonable values of the SR 
parameters of the order of 1 per cent. Let us note that the last equality 
in 0 is strictly valid for the modes very close to the pivot scale fc ~ = 

a^H. Away from the pivot scale small deviations proportional to the SR 
parameters, —2{esR — ? 7 s_R)ln^^j, are neglected. Depending on the SR 

parameters, on and on the vacuum choice {5k) the quantum corrections 
may lead to a power loss or a power increase for large scales which can be 
generically parametrized in the following form; 





(9) 


where the expression (|^ is simply obtained on comparing with ([^. An anal¬ 
ogous parametrization holds for the tensor sector with a different q. 


3.1 Extrapolation beyond NLO 

The parametrization of the primordial spectra by (|^ is still not suitable for 
comparison with observations. In the A; fc* limit the quantum gravitational 
corrections to po are either negative or very large (infinite in the /c —?• 0 limit). 
Such an apparently pathological behavior is simply the consequence of the 
perturbative technique employed to evaluate the corrections. One may hope 
that resummation to all orders leads to a finite result. In any case we are not 
allowed to extend the validity of the perturbative corrections up to O (1). 
Thus, instead of introducing a sharp cut-off on the NLO expressions for the 
modified spectra by multiplying q by an ad hoc step function which keeps 
the correction small but leads to a discontinuous spectrum, we interpolate 
our expression through a well defined function with a finite and reasonable 
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Table 1; Varying parameters 


In (lO^M,) 

Us 

r 

(Xg 

qi 

q2 

[2.7,4.0] 

[0.9,1.1] 

[0,0.8] 

[-0.1,0.!] 

|0,1| 

[0,0.5] 


behavior in the fc —?■ 0 limit. Such a function, which must reproduce ([^ 
when q {k^/kY 1, may be regarded as a resummation of the perturbative 
series. 

In order to restrict the number of parameters which will be htted by the 
comparison with the data and still allow for different limits when fc —>■ 0 we 
consider the following parametrization: 




(L)l + ?i (t)' 




(L) 

Po 


1 + (gi - q2) 



( 10 ) 


where one more parameter w.r.t. (|^ has been added in order to obtain a 
regular expression for k small. Let us note that the above modihcations are 
substantially different from considering a running spectral index as such as 





( 11 ) 


Indeed in the latter case the standard power law dependence is affected at 
both large and small scales and, in particular, a negative running would lead 
to a zero amplitude in the fc —)■ 0 limit and a smaller amplitude w.r.t. simple 
power law when k^ k^. The modihed spectrum (10) reduces to the power 
law case when k ^ k^ and may lead to a non zero amplitude when fc —)■ 0 
depending on the choice of the parameters qi^ 2 - 


4 Data Analysis 


In this section we report the comparison between the theoretical predictions 


given by (10) and Planck 2013 [Tj and 'Q\CEP2/Keck Array [H] datasets. 
The analysis is performed using the Markov Chain Monte Carlo (MCMC) 
code COSMOMC [13] which has been properly modihed to take into account 
the estimated quantum gravitational effects. 

Let us note that BD vacuum in the tensor sector (see |6] for more details) 
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gives a power increase at large scales in the tensor spectrum. Such an in¬ 
crease would be counterbalanced by a loss of power in the scalar sector as 
far as temperature correlations are concerned. One may parametrize such a 
power increase in a suitable way just as we did for the scalar sector in order 
to eliminate the divergence for small k and £t the corresponding parameter 
with the data at our disposal. Since our main source of data comes from 
temperature correlations, which do not discriminate between scalar and ten¬ 
sor fluctuations, we neglect a priori quantum gravitational corrections in the 
tensor spectrum. Such a choice is a simplifying assumption done in order not 
to have to disentangle possible degenerate parameters. Let us note, however, 
that such a choice can be realized physically either by an appropriate vacuum 
choice, differing from a pure BD, or by a very long cutoff scale associated 
with tensor dynamics. Thus we limit our analysis to a subset of the more 
general case for which the quantum gravitational corrections affect the ten¬ 
sor sector in a non negligible way. The tensor spectrum is then given by the 
unperturbed power law expression 


Pi = A 



( 12 ) 


and we assume that the LO spectra are generated by the conventional SR 
mechanism and single held inflation. As a consequence the consistency condi¬ 
tion relating scalar and tensor spectral indices and the tensor to scalar ration 
is valid. Indeed throughout the analysis we assume that the consistency rela¬ 
tion (already implemented in COSMOMC) between the spectral indices and 
the tensor to scalar ratio; 

«t = -^(2-«s-0 (13) 


holds to the second order in the SR approximation and consequently the 
amplitude of the spectrum of tensor perturbations is given by At = rAg. We 
then consider a primordial scalar spectrum p^^'> parametrized by 


Ps ^Po 




1+(T) 


(14) 


where 1 — 2^2 simply fixes the limit of Ps when A; —)■ 0. Let us note that ^2 = 0 
or gi = 0 correspond to the standard power-law case with no loss of power 
{Ps = Po^^) and q 2 = 0.5 corresponds to zero power at k = 0. The expression 
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Table 2: Models List 


Model 

Primordial spectra 

Datasets 

Parameters 

1 

2 

Power law 

no tensors 

PL 

PL+BK 


3 

4 

Power law 
and tensors 

PL 

PL+BK 

Ag, rig, r 

5 

6 

Running spectral index 
and tensors 

PL 

PL+BK 

Ag, rig, r, ag 

7 

8 

Quantum gravitational 
corrections and tensors 

PL 

PL+BK 

Ag, rig, r, qi, q 2 

is a parametrization equivalent to (10 

with qi = 

ql (1 - 2 ^ 2 ) and q 2 


which we have found to be more convenient to be used in COSMOMC. 

Our analysis is based on Planck+WP (PL) datasets released in 2013, which 
include the Planck TT data and the WMAP9 polarization data, and the 
BICEP2/AecA; Array dataset (BK) released in 2015. In particular we use 
CAMspec, lowLike, lowland /ensmg'publicly available Planck likelihoods. We 
End the best fit for our model with and without BK results and compare it 
with standard power law predictions and with predictions assuming a running 
spectral index ( [II| . In the standard power law case we either hx r = 0 or let 
it vary. 

For simplicity we obtained the best-fits for the parameters of the primordial 
spectra shown in Table and the parameters are taken to vary with uniform 
priors in the intervals indicated in the same table. The remaining parameters 
are fixed to the Planck best-fit in base_planck_ lowl_ lowLike_ highL run. Let 
us note that the pivot scale /c* is 0.05 Mpc“^ and is the same for both the 
scalar and the tensor sector. The priors for Ag, Ug, r and ag are those used 
by Planck 2013 analysis. The additional parameters qi and q 2 are chosen 
to vary in the largest possible interval reproducing a power loss for large 
scales (compared to the pivot scale) with the parametrization chosen. At 
present our theoretical predictions are not able to constraint the value of 
such parameters or estimate possible allowed intervals where to let them 
vary (see [I5] for an attempt to estimate priors from quantum gravity), thus 
the choice of broad enough priors seems reasonable. In particular the prior for 
qi is chosen in order to search for effects at scales larger that the pivot scale 
k~^, which seems reasonable given the present data, with the limiting value 
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Table 3: Monte Carlo Best-fits 


# 

In (101*^71,) 

Us 

r 

Q-s 

Qi 


1 

3.097 

0.957 

- 

- 

- 

- 

2 

3.097 

0.957 

- 

- 

- 

- 

3 

3.097 

0.958 

0.61 ■ 10-3 

- 

- 

- 

4 

3.097 

0.958 

0.24 • 10-3 

- 

- 

- 

5 

3.100 

0.957 

0.26 ■ 10-3 

-0.12 ■ 10-3 

- 

- 

6 

3.100 

0.957 

0.39 ■ 10-3 

-0.13 ■ 10-3 

- 

- 

7 

3.099 

0.954 

0.17 

- 

0.73 ■ 10-3 

0.13 

8 

3.098 

0.955 

0.44 ■ 10-3 

- 

0.59 ■ 10-3 

0.11 


gi = 0 which reproduces the standard power-law ps = paranietrization. 
The prior for q 2 is chosen to let the parameter vary between q 2 = 0 and 
q 2 = 1/2 where the standard power-law is obtained independently of qi. The 
values of q 2 with q 2 > 1/2, leading to an increase of power, and with q 2 < 0, 
leading to a physically unacceptable negative spectrum, are excluded from 
the analysis. 

The different combinations of primordial spectra and datasets considered 
are listed in Table and the corresponding fits are labelled with an index 
specifying the model number in table The best fits for the parameters we 
vary are presented in Table and the corresponding effective defined as 
—2 In £, where C is the likelihood, are listed in Tables|^and[^ In particular in 
the latter tables we present the effective derived from the full dataset used 
in each analysis {xTot) those derived from the comparison with only the 
Commander dataset or the BK dataset. Let us note that Commander data 
only regard temperature correlation (and not polarization) and are restricted 
to the lowest multipoles 2 < I < 50. Finally the differences between the total 
for the different cases are reported on using our model as reference. In 
particular the cases 1, 3 and 5 are compared with 7 and the cases 2, 4, 6 are 
compared with 8. 

4.1 Results 

The MCMC results (see Tables and show that the quantum gravita¬ 
tional modification of the standard power law form for the primordial scalar 
spectrum improves the fit to the data. Such improvements are much more sig- 
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Table 4: Monte Carlo Comparison (PL) 


# 

'X. C ommander 

XTot 

^Xi-X ^7 

1 

-5.56 

9864.2 

5.9 

3 

-5.61 

9864.2 

5.9 

5 

-8.72 

9860.6 

2.3 

7 

-10.6 

9858.3 

0 


Table 5: Monte Carlo Comparison (PL+BK) 


# 

X Commander 

x\k 

Xiot 

> 

bO 

III 

1 

00 to 

2 

-5.56 

46.9 

9911.3 

7.2 

4 

-4.66 

45.7 

9910.9 

6.8 

6 

-8.62 

45.3 

9906.7 

2.6 

8 

-10.6 

45.3 

9904.1 

0 



Figure 1: 68% and 95% confidence level constraints on r and Ug in cases 
5-6. 


nificant w.r.t the standard modifications of the primordial spectra obtained 
by considering a running spectral index. Concerning the best fit to the pa¬ 
rameters qi and q 2 we observe that the two datasets lead to close predictions 
with 

qi ~ 0.73 ■ 10"^ , 2q2 ~ 0.26 (15) 

when Planck data alone are considered and 

qi ~ 0.59 ■ 10“^ , 2q2 ~ 0.22 (16) 
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Figure 2: On the l.h.s. 68% and 95% confidence level constraints on r and 
Us in cases 7 — 8. On the r.h.s. 68% and 95% confidence level constraints on 
qi and q 2 in cases 7 — 8 . 


when BK data are added to the analysis. Correspondingly and Ag also 
take very similar values for the best £t. Let us note that a larger tensor to 
scalar ratio w.r.t. the case 5 with a running spectral index (see figure [I]j^ is 
allowed for case 7. 

The value of q 2 indicates a ~ 20 — 25% loss in power when k approaches 
zero and the small value of qi indicates that the quantum corrections are 
negligible around the pivot scale fc*. Let us note that the tensor to scalar 
ratio r may be quite large at the pivot scale when the analysis is restricted 
to PL data (but weakly constrained). The predictions concerning the value 
of r and the amplitude of the scalar perturbations Ag £x the energy scale of 
inflation 


Hi 

Mp 


— —Ag 
2 


(17) 


and in particular, on assuming r = 0.05 and In (lO^^As) = 3.1 (from the best 
£t of the case 8), one is led to ~ 0.05 ■ 10“®. The extrapolation (10) 

in the fc —)■ 0 limit, on comparing with ([^, gives the following constraint on 
the physical parameters of the model; 


fSk 17 


N4l-ns-r/S)\ 

18 ) 



(18) 


For reasonable values of iV* ~ 60, r 


0.05, qi and q 2 given by (15) and 


11 




















Figure 3: Marginalized 1-D likelihoods forr, qi and q 2 in cases 7 (black line) 
and 8 (red line). 


imposing the BD prescription for the vacuum, 6k = —k/2, we observe that 
the quantity inside the bracket on the hh.s. of (18) is ~ —2.5. Let us note 
that slightly different values for r, Us and iV* may lead to results which differ 
by an order of magnitude. 

Correspondingly, for the cases and ^8 we hnd the following estimate 
for the cut-off scale k 


k #7 ~ 38fc* ~ 1.9 Mpc ^, %8 — 30fc* ~ 1.6 Mpc ^ (19) 

which are very large compared to the wave number associated with the largest 
observable scale in the CMB namely fcmin ~ 1.4 ■ 10“"^ Mpc“^. Let us note 
that the existence of such a small fundamental scale may have relevant conse¬ 
quences on astrophysical observation. Indeed it is associated with distances 
which are comparable with the diameter of our galaxy. 

The primordial scalar spectrum with a running spectral index (case 6) 
is almost indistinguishable from case 8 in polarization (see Table where 
the effective BK data is presented) but is disfavored (Ay^ ~ 2.6 on 

adding only one free parameter to the £t) when Planck and BK data are 
simultaneously considered. In £gure|^we hnally plot the marginalized likeli¬ 
hoods for r, qi and g 2 - The correspondent marginalized 68%, 95% and 99% 
confidence intervals are listed in tables [6]j^ for cases 5, 6, 7 and 8 respectively. 


12 






Table 6; Marginalized confidence intervals - Case 5 



68% 

95% 

99% 

r 

[0,0.94-10-^] 

[0,0.18] 

[0, 0.24] 

— 10 ■ as 

o 

00 

o 

1 

o 

[0.22 ■ IO-SO.28] 

[-0.16 ■ 10-^0.32] 


Table 7: Marginalized confidence intervals - Case 6 




68% 

95% 

99% 

r 

[0.90 ■ 10 

0.69- 10-1] 

[0,0.10] 

[0,0.13] 

— 10 ■ as 

[0.67 ■ 

10-1,0.18] 

[0.94-10-2,0.25] 

[-0.23 ■ 10-1,0.28] 


Table 8: Marginalized confidence intervals - Case 7 



68% 

95% 

99% 

r 

[0, 0.24] 

[0, 0.47] 

[0,0.59] 

Qi 

[0.47-10-1,0.10] 

[0.66 - 10-2,0.15] 

[0, 0.24] 

0.2 

[0.63 - 10-1,0.19] 

[0.45 - 10-2,0.26] 

[0,0.30] 


Table 9: Marginalized confidence intervals - Case 8 



68% 

95% 

99% 

r 

[0.12-10-1,0.79-10-1 

[0,0.11] 

[0,0.14] 

Oi 

[0.23 - 10-1,0.10] 

[0,0.19] 

[0,0.39] 

02 

[0.35 - 10-1,0.13] 

[0,0.18] 

[0, 0.24] 


The marginalized likelihoods for qi and q 2 show a 2a deviation from 
standard power law for case 7 and a la deviation for case 8. Let us note that 
the marginalized likelihood for q 2 rules out q 2 = 1/2 (i.e. the case of no /c/fc* 
dependence at the numerator in (14)) at more than 3a in both case 7 and 8. 
Thus the comparison with the data leads to severe restrictions on the form 
of the quantum gravitational corrections (10) as a function of {k/k)^. 


5 Conclusions 

The possibility of a loss of power at large scales in the spectrum of the CMB 
temperature anisotropies has intriguing consequences. On the phenomeno- 
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logical side the coexistence of these effects signals the necessity of modifying 
the standard power law assnmption for the primordial spectra predicted in 
the context of SR inflation. On the theoretical side both these effects can 
be related to the physics at and above the Planck scale and (even if still 
nncertain given the present precision of the data) shonld be interpreted in 
the context of quantum gravity. 

In this paper we present an analysis of the effects of the quantum gravita¬ 
tional corrections on the spectra of primordial perturbations calculated in the 
context of a Wheeler-De Witt approach. Such effects have been obtained and 
discussed in detail in [7] and have a quite distinctive form which affects more 
the large scales of the CMB. For simplicity the analysis has been restricted 
to the particular case of negligible quantum gravitational contributions on 
the spectrum of primordial gravitational waves. The NLO predictions for 
the scalar sector have been described in terms of two parameters, suitably 
extrapolated beyond the NLO and examined down to fc —)■ 0. Finally, we 
compared PL and BK observation with our predictions by using a common 
parametrization for the primordial spectra and varying only a few relevant 
parameters in the cases discussed. Other parameterizations have been con¬ 
sidered, however the one we present is the simplest and leads to the best 
results. 

The results show that the effect of our quantum gravitational corrections 
compared to the standard cases leads to a better £t with the data. A similar 
result has been obtained in the literature |5] by adopting a different modih- 
cation of the primordial spectra leading to a power loss at large scales and 
a similar improvement to the £t was found. Nonetheless the possibility of 
explaining the loss in power as an effect of quantum gravity on standard 
SR inflation is intriguing and the results obtained suggest the possibility of 
looking for quantum gravity effects at much smaller scales. On including 
the BK data in our analysis we found that the statistical relevance of the 
possibility of a loss in power for large scales is nearly invariant. Furthermore 
the comparison with the data predicts, for our model, a loss in power of 
about 20 — 25% w.r.t. the standard power law and hxes the scale k which 
necessarily appears in the theoretical model. If such a scale plays a crucial 
role in the description of gravitational phenomena one may expect that other 
astrophysical observables are affected around it. Unfortunately, due to the 
large error bars in the temperature correlations at large scales, observations 
are not able to severely constrain the parameter of the model. A la — 2a 
deviation from the standard power law parametrization is found. Still the 
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quite distinctive form of the quantum gravitational corrections gives a clear 
indication on their shape. 
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